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ABSTRACT 

We present Early Science observations with the Large Millimeter Telescope, AzTEC 1.1 mm 
continuum images and wide bandwidth spectra (73 — 111 GHz) acquired with the Redshift 
Search Receiver, towards four bright lensed submillimetre galaxies identified through the 
Herschel Tensing Survey-snapshot and the Submillimetre Common-User Bolometer Array-2 
Cluster Snapshot Survey. This pilot project studies the star formation history and the physi¬ 
cal properties of the molecular gas and dust content of the highest redshift galaxies identified 
through the benehts of gravitational magnihcation. We robustly detect dust continuum emis¬ 
sion for the full sample and CO emission lines for three of the targets. We find that one 
source shows spectroscopic multiplicity and is a blend of three galaxies at different redshifts 
{z = 2.040, 3.252 and 4.680), reminiscent of previous high-resolution imaging follow-up of 
unlensed submillimetre galaxies, but with a completely different search method, that conhrm 
recent theoretical predictions of physically unassociated blended galaxies. Identifying the de¬ 
tected lines as ^^CO (J„p = 2 — 5) we derive spectroscopic redshifts, molecular gas masses, 
and dust masses from the continuum emission. The mean H 2 gas mass of the full sample 
is (2.0 ± 0.2) X IO^^Mq/p, and the mean dust mass is (2.0 ± 0.2) x Mq/ ft, where 
/r w 2 — 5 is the expected lens amplihcation. Using these independent estimations we infer a 
gas-to-dust ratio of Jgdr ~ 55 — 75, in agreement with other measurements of submillimetre 
galaxies. Our magnihed high-luminosity galaxies fall on the same locus as other high-redshift 
submillimetre galaxies, extending the ■^co -Lfir correlation observed for local luminous and 
ultraluminous infrared galaxies to higher far-infrared and CO luminosities. 

Key words: submillimetre: galaxies - galaxies: high redshift - galaxies: evolution - galaxies: 
ISM - cosmology: observations 
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1 INTRODUCTION 

Observations of nearby galaxies with the Infrared Astronom- 
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(FIR)/submillimetre background with the Cosmic Background Ex¬ 
plorer (COBE), showed that the Universe emits a comparable en¬ 
ergy density at infrared (IR) and submillimetr e wavelengths as it 
does at optical and ultraviolet waveb ands fe.g. ISoifer et al.iri987l : 
IPuget et ^Il996l : Ipixsen et aDll998h . A breakthrough in resolv¬ 
ing this background came with the discovery of a large popula¬ 
tion of bright sources at high redshift through sin gle-dish telescopes 
observations at submillimetre wavelengt hs le.g Ismail et ‘ill 19971: 
ISarger et"^ 1 19981 : iHughes et alj Il998h . These (sub)millimetre- 
selected galaxies (hereafter SMGs) are characterized by large FIR 
luminosities (> 10^^ Lq), large star formation rates (SFRs, > 300 
Mq yr“^), large gas reservoirs (> 10^° Mq), and a number 
density that is high co mpared to local u ltra-luminous IR g alaxies 
dChapman et alj2005l : see also review bv ICasev et al.ll20i^ . 

However, identifying and understanding the nature of these 
sources has proven to be challenging due to the low angular reso¬ 
lution of single-dish telescopes and the faintness of these galax¬ 
ies in the rest-frame optical and ultraviolet bands, w hich makes 
the a ssociation with the correct counterpart difficult dSlain et ^ 
l2002h . High-resolution radio/mid-infrared observations have often 
been used to identify counterpart galaxie s on 1 — 2 arcsec scales 
(e.g. Ilvison et aHl2002l : IPope et aPbOOSh . providing the precision 
needed to target optical redshifts, which are essential to understand 
the physical nature of these sources. However, recent submillimetre 
interferometric observations have shown that resourcing to surro¬ 
gate wavelengths could miss ~ 45 per cent of SMGs, and of those 
which are claimed to be iden t ified, approx i matel y one third are 
incorrect dSmolcic et alJ 1201 2l : iHodge et al.l l2013l : ISimpson et al.l 
l2015t) . Furthermore, the use of these wavelengths suffers from a 
well-known systematic bias against high-redshift (« > 3 ) sources 
and therefore our understanding of these galaxies could be incom¬ 
plete. 

Millimetric spectroscopy provides the best way to determine 
redshifts for this population of galaxies, since they are expected to 
have luminous emission lines in the millimetre bands. Furthermore, 
the cold molecular phase of the interstellar medium (ISM) is an ex¬ 
cellent tracer of the molecular gas reservoir, and the overall galaxy 
dynamics. For instance, CO(l-O) has been extensively used to map 
the molecular gas conten t of this kind of gal a xies (e.g.|CariIli_eLal,| 


20ld: [Harris et alj|20ld: [i vison et ^l20 ld: iRiechers et al. 
Bothwell et alj 2013 ; see also review by Carilli & Waltei 
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120131). 


When combined with the SFRs (or FIR luminosities), gas masses 
also provide estimations of the star formation efficiencies (SFE), 
whi ch is an important parameter for models of galaxy evolution 
(e.g. lSantini et al.ll20l4lDessauges-Zavadskv et alj|2015h . 

However, CO spectroscopy is very time consuming, when re¬ 
tuning receivers to search for molecular line emission, requiring 
entire nights to detect a si ngle source, and in some cases wit h 
no successful detection (e.g. iGreve et^l2005l : iHuvnh et al.ll2014h . 
This method has only recently become competitive with the in¬ 
creased bandwidth of (sub)mm facilities, but is still not yet suffi¬ 
ciently efficient to obtain redshifts for su bstantially large samples 
of unlensed SMGs jBothwell et al.ll2oT^ . On the other hand, CO 
line observations for strongly lensed systems can be obtained more 
easily, represe nting one route to study the molecular gas content 
of S MGs (e.g. Knudsen et alj|200^ : ICox et ai]l201ll : [Harris et"^ 
[2012l:[WeiBetal.[[2013h . Although in some cases, differential mag¬ 
nification could affect the pro perties derived from observations 
(e.g. line ratios. l^rieantllloTil) . cluster lensing provide us more 
reliable results since is rare to find this kind of magnification gra¬ 
dients. 

Here we present new observations of a sample of four cluster- 


lensed SMGs, discovered by the Herschel Lensing Sur vey-snapshot 
(HLS -snapshot, Egami et al. in preparation; see also [Egami et al.[ 
[mi3) and subsequently detected with the Submillimetre Co mmon- 
User Bolometer Array-2 (SCUBA-2, [Holland et alj[2013l) on the 
James Clerk Maxwell Telescope (JCMT) as part of the SCUBA- 
2 Cluster Snapshot Survey (S2CSS, Geach et al. in preparation). 
The new observations comprise imaging from the 1.1 mm contin¬ 
uum camera AzTEC jWilson et alj[200^ and the wide bandwidth 
(73 — 111 GHz) sp ectrometer Redshift Search Receiver (RSR, 
[Erickson et alj|200% ins talled on the Large Millimeter Telescope 
Alfonso Serrano (LMlBlHugher ^t_^lJl20J^3) > located on the sum¬ 
mit of Volcan Sierra Negra (Tliltepetl), Mexico, at an altitude of 
~ 4600 m. Combining the high sensitivity and bandwidth of the 
LMT instruments with the higher flux density of amplified galax¬ 
ies, the detection of CO molecular lines and dust continuum can 
be reached in reasonable times, allowing us to derive spectroscopic 
redshifts, and therefore to inquire into the dust and gas properties 
in these galaxies. 

In §2, we describe the sample selection, observations and data 
reduction. In §3, we describe the line identification, redshift estima¬ 
tions, and gas and dust properties. In §4, we discuss the implication 
of these results for the gas-to-dust ratio and SFE. Our results are 
summarized in §5. 

All calculations assume a A cold dark matter cosmology 
with Oa = 0.68, Qm = 0.32, and Ho = 67 kms“^Mpc~^ 
jPlanck Collaboration et alj2013l) . 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample selection and SCUBA-2/JCMT data 


The four bright SMGs targeted here were originally discovered 
by the HLS-snapshot survey, which has obtained shallow but 
nearly confusion-limited Spectral and Photometric Imaging Re¬ 
ceiver (SPIRE) images of 279 Massive Cluster Survey (MACS) 
clusters. Three of the four galaxies have a SPIRE flux densities 
of more than 100 mjy at either 250, 350 or 500 ^m, and have sub¬ 
sequently produced bright detections in the S2CSS. The remaining 
target (HLS J102225.9-1-500536) was detected by HLS below the 
100 mJy threshold but was comparably bright at 850 p.m. The anal¬ 
ysis of the Herschel data and their optical-IR counterparts will be 
presented elsewhere by the HLS team. 

S2CSS exploited poor weather conditions at the JCMT (opac¬ 
ities of r 225 GHz = 0.12-0.2) to search for bright (Sgsoum > 30 
mJy) gravitationally lensed S MGs, comparable to the ‘Cosmic Eye¬ 
lash’ Iswinbank et al.[|2010[) . in the cores of hundreds of rich clus¬ 


ters of galaxi es, primarily selecte d from the MACS and SMACS 
survey s (e.g. Ebeling et al.| [200l[ : see also target selection of the 
HLS in tegami et al.l2010h . The ‘Cosmic Eyelash’ is a typical SMG 
at 2 « 2.3 with an intrinsic flux of Ss 7 o ~ 3 mJy, which is magni¬ 
fied 30 times by a foreground cluster, appearing as an Sgro ~ 106 
mJy source. The high magnification and spatially-resolved struc¬ 
ture allows a very detailed study of the internal structure and chem¬ 
istry of this galaxy on < 100 pc scales. Although rare, the bright 
nature of these sources makes them detectable at high significance 
even when observing conditions are relatively poor. Eor each clus¬ 
ter target a 30 minute DAISY scan map was obtained with SCUBA- 
2, which allowed the detection of luminous SMGs within a few 
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Figure 1. From left to right, 30 X 30 arcsec^ Infrared Array Camera (IRACySpteer images in the 3.6 fim band at the SCUBA-2 position of HLS 
J045518.0-H070103 and HLS J204314.2-214439 andHSr images in the F606VF band at the position of HLS J102225.9-H500536 and HLS J173140.8-H225040. 
The SCUBA-2 850 fim and AzTEC 1.1mm S/N contours are also shown in green (solid line) and red (dashed line), respectively. Contour levels stall at 1.5 (t 
and are spaced in steps of 1 .5 it. See the elech'onic edition of the journal for a colour version of this figure. 


Table 1. LMT observing log for data taken during the Early Science Phase in 2014. Column 1, name of the targeted galaxies; columns 2, 3, and 4: dates, 
opacities and integration times of AzTEC observations; columns 5, 6, and 7: dates, opacities and integration times of RSR observations. 


ID 

Date 

AzTEC 

Opacities 

(t225GHz) 

tint 

(min) 

Date 

RSR 

Opacities 

('r22SGHz) 

tint 

(min) 

HLS J045518.0-t070103 

Feb 10 

0.12 

20 

Feb 19, Mar 02, Apr 02 

0.06, 0.17, 0.26-0.29 

45 , 45,20 

HLS J102225.9-t500536 

Jan 31 

O.IO 

20 

Mar 27, 29, Apr 02, 07 

0.17, 0.15,0.20, 0.15 

10 , 60, 50, 35 

HLS JI73140.8-t225040 

Mar 21, May 18 

0.12, 0.18 

4,10 

Apr 21, May 03 

0.30, 0.30 

30,75 

HLS J204314.2-214439 

May 14, 16 

0.20, 0.22 

10,10 

May 08, 12,18, 25,26 

0.46-0-57, 0.29-0.32, 
0.13-0.15,0.42-0.44, 0.35 

90, 30, 50, 60, 30 


arc minutes from the cluster cores (depending on source bright¬ 
ness). The data were r educed with the Sub -Millimetre User Reduc¬ 
tion Facility (SMURF; Chapin e t al .11 20131) , following a procedure 
described in lOeach et al.l ll2013l) where the map making was opti¬ 
mized for point-source detection. Typical root mean squared noise 
(rms) noise in individual 850-jrm maps was 5-10 mjy (see Table 

H. 


We conducted observations with AzTEC and RSR on the LMT 
towards a sample of four 850 fim sources that represent some 
of the most secure (> 5 (t) detections within S2CSS around the 
clusters MACS J0455.2-I-0657, Abell 980 (MACS J1022.4-b5006), 
MACS J1731.6-b2252, and RXCJ2043.2-2144 (MACS J2043.2- 
2144) at g = 0 .425, 0.158, 0.366, and 0. 204, respectively 
tebeling et alJl99i.|201(]|; iBdhringer et al.ll200l; iMann & Ebelind 

I2OI2I) . Postage stamps of Spitzer or Hubble Space Telescope 
(HST) images, and overlaid with SCUBA-2 and AzTEC contours 
at these positions are shown in Fig. [T] Details of these bright 
galaxies, HLS J045518.0-b070103, HLS J102225.9-I-500536, HLS 
J173140.8-1-225040, and HLS J204314.2-214439, and their fluxes 
are given in Table 


tions were conducted over several observing nights with an opacity 
range of r 225 GHz = 0.07 — 0.18 and a typical integration time of 
~ 20 min per source, achieving an rms of 2-3 mJy (see Table[T]( in 
the individual maps. 

Each individual AzTEC observation is composed of a set of 
timestreams which store all the bolometer sky signals, the projected 
position of the array on the sky and the weather conditions. Cali¬ 
bration is performed using observations of the asteroids Ceres and 
Pallas. To generate the final maps from the ra w data we use th e 
AzTEC Standard Pipeline described in detail bv IScott et alj ( l2008h . 

In Table we summarize the 1.1 mm photometry measure¬ 
ments from the AzTEC maps shown in Eig. The AzTEC po¬ 
sitions, which, for these observations, have a typical uncertainty 
of 2 arcsec, are consistent with the SCUBA-2 positions within 
the 95 per cent confidence intervals, except for the case of HLS 
J204314.2-214439. However, as discussed in 113.1.41 our RSR 
spectrum has revealed three sources in the beam for this target and, 
therefore, the estimated positional uncertainty could be underesti¬ 
mated. 


2.2 AzTEC/LMT 

Observations were obtained using the 1.1 mm continuum camera 
AzTEC on the LMT between 2014 January and May (see Table 
[T](. During this Early Science Phase operation, only the inner 32-m 
diameter section of the telescope surface is illuminated, leading to 
an effective beam size of 8.5 arcsec, a factor of 2 better than 
the SCUBA-2 850 /rm-beam. The scanning technique provided a 
Lissajous pattern covering a 1.5 arcmin diameter region with uni¬ 
form noise in the map centred at the SCUBA-2 position. Observa- 


2.3 Redshift Search Receiver/LMT 

Observations with RSR were obtained between 2014 February and 
May. The RSR has 4 pixels arranged in a dual beam, dual polar¬ 
ization configuration. The four broad-band receivers cover instan¬ 
taneously the frequency range 73 — 111 GHz with a spectral res¬ 
olution of 31 MHz (~ 100 km s“^ at 93 GHz). The current 32-m 
telescope gives us an effective beam size of 20 arcsec at 110 GHz 
and 28 arcsec at 75 GHz. 

Observations were obtained at the SCUBA-2 positions of the 
four galaxies listed in Table|^ over several observing nights with an 
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Table 2. Summary of SCUBA-2 and AzTEC photometry. Colum n 1: name of the tar geted galaxies; columns 2, 3, 4, and 5: position, 95 per cent confidence 
interval for positional uncertainties following the description by llvison et for single sources, signal-to-noise ratio, and 850 /im flux density for 

SCUBA-2 obseiwations; columns 6, 7, 8, and 9: same information for 1.1 mm AzTEC observations. 


ID 

SCUBA-2 position 
(J2000) 

A(a,5) 

(arcsec) 

S/N 

SsSO/im 

(mJy) 

AzTEC position 
(J2000) 

A{q,(5) 

(arcsec) 

S/N 

Si.1mm 

(mJy) 

HLS J045518.0+070103 

04:55:18.1, +07:01:01 

3.2 

7.0 

37 + 5 

04:55:18.1, +07:01:05 

1.6 

8.0 

23 + 3 

HLS 1102225.9+500536 

10:22:26.6, +50:05:39 

4.1 

5.5 

45 + 8 

10:22:26.2, +50:05:41 

2.2 

5.9 

11 + 2 

HLS 1173140.8+225040 

17:31:40.5, +22:50:35 

2.6 

8.5 

49 + 6 

17:31:40.6, +22:50:37 

1.3 

9.8 

25 + 3 

HLS 1204314.2-214439 

20:43:14.5, -21:44:37 

1.7 

13.0 

66 + 5 

20:43:14.2, -21:44:43 

1.5 

8.4 

29 + 3 


AzTEC 



-3 0 3 5 8 

S/N 


RSR HLS J045518.0 + 070103 



HLS J102225.9 + 500536 



HLS J173140.8 + 225040 



HLS J204314.2-214439 



Figure 2. Left: 40 x 40 arcsec^ AzTEC/LMT 1.1 mm S/N maps. All the maps are centred at the AzTEC positions listed in Table[2| and have been scaled for 
display purposes. Right: RSR/LMT spectra of the four HLS/S2CSS targets observed with the LMT. The RSR/LMT spectra have been rebinned into 3 pixels 
bins for better visualization. We mark with horizontal lines the 2 and 3.5(7 average noise level for each spectrum, where 1(7 has been calculated as the standard 
deviation of the whole binned spectrum after removing the identified lines. The transitions detected above the adopted threshold detection are marked with 
vertical dashed lines. In the spectrum of HLS J204314.2-214439 we have identified contributions from three different galaxies at different redshifts. We also 
mark with an arrow the expected position for the [Cl] transition at 2 = 4.680 (see ^3.1.4L See the electronic edition of the journal for a colour version of this 
figure. 
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opacity range of r 225 GHz = 0.13 — 0.57 (see Table[T) and system 
temperatures of Tays = 90 — 140 K with an average Taya ~ 100 K. 
Pointing was performed every hour on bright millimetre sources, 
leading to a pointing accuracy (rms) ~ 2 arcsec. 

The data were reduced and calibrated using Data REduction 
and Analysis Methods in PYthon (DREAMPY), which is the RSR 
data reduction pipeline software written by G. Narayanan. After 
flagging any data adversely affected by a hardware or software 
problem, a linear baseline is removed from each spectrum. The final 
spectra were obtained by averaging all scans using 1 /(t^ weights. 
The spectra in antenna temperature units (T^) are converted to flux 
density units using the conversion factors 6.4 and 7.6 Jy K“^, for 
73 < Si 92 GHz and 92 < jz < 111 GHz, respectively, which 
are based on the calibration of Uranus and MWC349A, observed 
at EL = 35° — 75°, conducted in the Early Science Phase. Pinal 
spectra are shown in Fig.[2 

2.4 SMA observations 

Observations on the cluster RXCJ2043.2-2144 were obtained with 
the Submillimeter Array (SMA) as a follow-up programme of the 
Herschel Lensing Survey (HLS) on 2012 May 8. After the LMT 
obtained a secure redshift for a z = 4.68 galaxy in this field (see 
113.1.4t . these data were reanalysed to search for the predicted [CII] 
line emission at t/ = 334.6 GHz as a means of further confirmation. 

Observations were made using six antennas in a compact con¬ 
figuration. An on-source integration time of 3.9 h was obtained with 
the primary LO set to 342.046 GHz, using a single polarization SIS 
junction receiver with 4 GHz bandwidth per sideband. This band¬ 
width covered the expected frequency of the redshifted [CII] line 
for z = 4.68 (rest frequency 1900.54 GHz shifted to 334.60 GHz). 
The weather was very good, with the atmospheric opacity at 225 
GHz varying from 0.05 to 0.07 throughout the observations. The 
synthesized beam size was 2.4 x 1.9 arcsec^. The primary flux 
density calibrator was Neptune, which provided a flux scale accu¬ 
rate to ~ 5%. Instrumental and atmospheric gains were calibrated 
using complex gain calibrators J1924—292 (4.05 Jy), J2000-178 
(0.46 Jy), and J2158-150 (1.66 Jy). The visibility data were resam¬ 
pled and averaged to a velocity resolution of 20 km s“^, and then 
binned further to 100 km s“^ resolution in the imaging routine. The 
rms noise obtained in the final images was 15.2 mjy beam“^ per 
100 km s~^ channel. 

These observations showed the presence of three distinct 
sources with a clear emission line feature at 334.6 GHz. Averag¬ 
ing consistent the spectra from these three positions produced the 
spectrum shown in Fig. |5] (rms of 8.8 mJy), which we identify as 
the [CII] transition. The analysis of the triple system and constraints 
to the lensing model of RXCJ2043.2-2144 (MACS J2043.2-2144) 
will be presented in a subsequent paper. 


3 RESULTS 

3.1 RSR analysis: line identification and spectroscopic 
redshifts 

Since the main goal of these observations is to blindly search for 
molecular lines in our spectra, with no previous knowledge on the 
expected observed frequencies for these lines, it is important to de¬ 
fine a criteria S/N threshold (^thresh) that ensures the reliability of 
detections and that maintains the false detection rate very low. 

To estimate the global noise of the spectra we fit the histogram 



- 1.0 - 0.5 0.0 0.5 1.0 

T • [mK] 

Figure 3. Histogram of pixel values in the RSR/LMT spectrum of HLS 
J102225.9+500536, as an example of the noise properties in the spectra. 
The histogram is well fitted by a Gaussian function, which is represented 
by the solid line. These Gaussian functions are used to make mock noise 
spectra in order to define a reliable threshold line detection. 


of pixel values with Gaussian distributions. As can be seen in Fig. 
[3 a Gaussian fit reproduces very well the spectrum, and can be 
used as an estimation of the average noise. This is true for all the 
acquired spectra, achieving fits of x?ed = 0-57 — 0.96. Using these 
Gaussian distributions, we make mock noise spectra in order to de¬ 
fine a robust detection threshold. These simulations assume that the 
noise is uniform at all frequencies and uncorrelated between pix¬ 
els. Using 100,000 realizations we find that the probability to have 
two adjacent pixels with S/N^ 3.5 is less than 0.1 per cent, and 
therefore this S/N ratio could be used as a detection threshold. To 
investigate the reliability of these simulations we search for neg¬ 
ative peaks, which are expected to be just noise, that satisfy this 
S/N threshold. We find none, either in the raw or in binned spec¬ 
tra, which supports the adoption of ^thresh = 3.5 (in two adjacent 
pixels) as a robust detection threshold. 

In order to identify line candidates in our spectra, we search 
for peaks with S/N^ ^thresh, where the noise is estimated as the 
standard deviation of the whole spectrum. If a peak is identified, we 
mask a region of 10 pixels centred at the position of the peak (since 
we know that the full width at half-maximum (FW HM) of lines de¬ 
tected in SMGs can be up to ~ 1000 km s"*^, e.g. lSothwell et al.l 
l2013h . This process is repeated until there are no more pixels with 
S/N^ ^thresh. If a line candidate is identified, we inspect it visually 
to confirm that it does not correspond to a bad baseline. In all our 
spectra we do not find any such case. In order to find weaker emis¬ 
sion lines, we rebin our spectra into 200 km s”*^ pixels (2 pixels 
bins) and repeat the process again. Once we have the position of 
line candidates we carry out a Gaussian line profile fit in order to 
obtain the central frequency, peak flux, and line width. 

Although the general noise is well fitted by a Gaussian distri¬ 
bution (see Fig. [3 and the baseline fluctuations in each spectrum 
are well confined within the ±2cr noise range (see Fig. |^, there 
are some baseline residuals that exceed this range or where the 
noise seems correlated. These baseline features are usually broad 
(> 1000 km s”*^) and therefore could preclude the detection of 
low signal-to-noise ratio (S/N< 5) broad lines (albeit the fraction 
of SMGs with very b road lines, > 1000 km s”*^, have been found 
to be low, e.g. 1/32. ISothwell et ai]|2013h . Broad lines could be 
missed since our method is based on the detection of high S/N 
pixels. However, a different method has been developed, which ex¬ 
ploits the full spectral information present in the RSR data by cross- 
correlating the observed spectrum with a theoretical or an empiri- 
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Table 3. Observational parameters of the lines detected in the RSR spectra, using a blind search. The three blended sources in the HLS J204314.2—214439 
target are identified with the extensions /SI, /S2, and /S3. Column 1: name of the source; column 2 and 3: signal-to-noise ratio; column 4; central frequency; 
column 5; identified transition; column 6; redshift; column 7; integrated flux; column 8; FWHM of the line; column 9; notes. 


Source 

S/N“ 

S/N*’ 

^obs 

(GHz) 

Transition 

z 

ScoAU 
(Jy km s“^) 

FWHM 

(km 

Notes 

HLS J045518.0+070103 

4.5 

7.3 

88.051 

i2C0(3-2) 

2.927“ 

4.7 ±0.8 

480 ± 70 

Alternative solution 2 : = 1.618 

HLS J173140.8+225040 

4.1 

5.7 

93.163 

12cO(3-2) 

2.712 

2.3 ± 0.8 

420 ± 120 

2 confirmed by IRAM 

HLS J204314.2-214439/S1 

10.3 

15.5 

75.827 

12 cO(2 -1) 

2.040 

6.1 ± 0.9 

530 ± 60 

2 confirmed by IRAM and VLT 

HLS J204314.2-214439/S2 

2 .1'* 

2 .0^* 

81.328 

12cO(3-2) 

3.252 

0.5 ± 0.3 

260 ± 130 

2 confirmed by VLT 


5.1 

5.7 

108.428 

12cO(4-3) 

3.252 

1.3 ± 0.4 

190 ± 40 


HLS J204314.2-214439/S3 

7.7 

12.3 

81.178 

12cO(4-3) 

4.680 

4.8 ± 0.6 

570 ± 60 

2 confirmed by SMA 


5.1 

8.0 

101.459 

12cO(5^) 

4.681 

2.8 ±0.5 

410 ± 50 



“The maximum S/N of the peak pixel from either the raw or the 200 km s“^ binned spectra. In this case, the noise has been estimated as the standard deviation 
of the whole spectnim after removing the identified lines. *’The S/N of the integrated line in the raw spectra (100 km s“^ pixels). “Tentative solution. “^This line 
is not formally detected by our algorithm, but searched for after the VLT redshift solution was proposed. 


cal spectral template (see lYun et al]|2007h . A detailed description 
of this method an d its applic ation to the SMG COSMOS AzTEC- 
1 is described by lYun et all (1201 Sl) . and we briefly summarize it 
here. A cross-correlation product is derived as a function of red¬ 
shift from the observed spectrum and a model spectral template. 
This model is weighted by a function which represents the relative 
strength of different molecular transitions, and for this, an empir¬ 
ical composite spectrum base d on observed relat ive line strengths 
for high redshift sources (e.g.. lSr)ilker et alj|2014h is adopted. This 
is a powerful method to derive redshifts which takes into account 
the information of many transitions which are not individually de¬ 
tected with a good S/N. We find that the redshift results are robust, 
and both methods give consistent results, when the integrated S/N 
of the lines are sufficiently high (S/N> 5). 

In Table[3]we summarize the lines detected in our spectra and 
Fig. a shows them individually. All these detections have at least 
one peak pixel at S/N^ 4.0 in the raw or in the binned spectra and 
at least one adjacent pixel with S/N^ 3.5 (except for a line that 
is not formally detected by our algorithm, see l]3.1.4b . However, a 
better estimation of the reliability of these detections is the S/N of 
the whole integrated line, which is also reported in Tabled All the 
lines detected originally by our algorithm have integrated S/N> 5. 

In order to derive the spectroscopic redshifts we have first to 
associate the detected lines candidates with the correct molecu¬ 
lar transitions. The simultaneous frequency coverage of the RSR 
means at least one CO transition falls within the spectral cover¬ 
age at all redshifts except for a narrow redshift range (the ‘red¬ 
shift desert’) in the range 0.58 < a < 1.08. Furthermore, two or 
more molecular line transitions fall within the RSR spectral range 
at « > 3.15 if we consider the possibility of a [Cl] line detection 
in addition to the CO lines. On the other hand, since we know that 
the spectral line en ergy distribution (SLED) peaks at Jup = 6 — 5 
for typical SMGs dSothwell et al.l 1^131 ; ISnilker et al.ll2014h . the 
absence of lines could be used to discard the high-redshift solu¬ 
tions for which we expect more than one line in our spectra. More¬ 
over, our selection of bright sources through galaxy clusters, im¬ 
plies that these very bright SMGs are expected to be gravitation¬ 
ally lensed. This selection tends to exclude low-redshift solutions, 
a < 1 — 1.5, due to the low probab i lity of being lensed at these 
redshifts, as discussed by lWeiB et alj l l2013l) . Taking these consid¬ 
erations into account, we associate our line candidates with differ¬ 
ent (Jup = 2 — 5) transitions and derive the spectroscopic 


redshift for each galaxy (see Tabled- Below we describe briefly 
the analysis for individual targets. 


3.1.1 HLS J045518.0+070103 


A line with an integrated S/N=7.3 with a well-fitted Gaussian pro¬ 
file (see Fig. 0 has been detected in this spectrum. The detection of 
Just one line al u = 88.051 GHz gives a tentative redshift for this 
galaxy and therefore we have to associate this line with the most 
probable transition based on other characteristics of the galaxy. If 
we associate this detection with CO(4-3) at 2 = 4.236, we ex¬ 
pect to detect the CO(5^) transition at t/ ~ 110.1 GHz still within 
our bandpass. The spectrum, however, does not show any line fea¬ 
ture around this frequency and, therefore, this solution is discarded. 
Other possible solutions are CO(3-2) or CO(2-l), which corre¬ 
spond to redshifts of 2 = 2.927 and 2 = 1.618, respectively. As we 
expect that th e probability for these galaxies to lie around 2 = 2 — 3 
is higher (e.g. Chapman et al.l200^ ; lYun et al.l201 j ; ISimpson et al.l 


l20li ; IZavala et alj|2014l) . we have associated our line candidate 

with the CO(3-2) at 2 = 2.927, although we cannot rule out the 
alternative lower redshift solution. 


3.1.2 HLS J102225.9+500536 

For this source we have no detections at S/N)5 3.5, either in the 
raw or in the binned spectrum. One possible explanation is that the 
galaxy lies in the ‘redshift desert’ of the RSR. However, as we can 
see in Table 0 this galaxy has the lowest flux density at 1.1 mm, 
which is less than half of those of the rest of the sample. This would 
suggest that we need to integrate more in this fainter galaxy to de¬ 
tect a CO emission line. Furthermore the ratio S'ssoMm/S'i.imm is 
the largest of the sample, which could be due to a higher tempera¬ 
ture, and make more difficult the detection of molecular lines than 
in the other sources. The upper limit for a line could be calculated 
using the equation Ico = 3(T(5nAnFWHM)^^^, where cr is the rms 
noise, Sv the velo city resolution, and Awfwhm the line width (e.g. 
iGreve et^l2005h . We adopted a line width of 500 km s“^ similar 
to the lines detected in the remaining spectra (see Table0, and a a 
derived from the standard deviation of the whole spectrum. Using 
this equation the upper limit for a line is 1.2 Jy km s~^. 
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HLS J045518.0 + 070103 


HLS J173140.8 + 225040 




HLS J204314.2-214439/S1 


HLS J20431 4.2-21 4439/S2 + /S3 




HLS J204314,2-214439/S2 


HLS J20431 4.2-21 4439/S3 



108,0 108.2 


108.4 108.6 

Frequency [GHz] 



HLS J204514.2-214439/S3 



Figure 4. LMT raw spectra at the position of the and [Cl] lines detected in our RSR spectra, best-fitting Gaussian profiles, and redshift derived from 

each individual line. The three blended sources in the HLS J204314.2-214439 target are identified with the extensions /SI, /S2 and /S3. The [CI](^Pi Pq) 
transition at ^ = 4.680 for HLS J204314.2-214439/S3 is only detected at Scr, and the FWHM of the Gaussian fit is consistent with the line widths of the two 
transitions. 
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3.1.3 HLSJ173140.8+225040 

This is another case of single-line detection at i/ = 93.163 GHz 
with a line integrated S/N=5.7. In this case we can associate this 
transition with CO(4-3) at « = 3.949 or CO(3-2) at « = 2.712 
having rejected the CO(2-l) solution at « = 1.475 due to the 
low probability of lensing at this lower redshift. Independently, 
earlier observations made by the HLS team had already identified 
the redshift of this galaxy through a CO redshift search with the 
IRAM 30-m/EMIR, yielding « = 2.712 through the detections of 
CO(5^) and CO(3-2) (Dessauges-Zavadsky et al., in preparation), 
and hence our single line detection unambiguously corresponds to 
CO(3-2). 

The line-shape shows some departure from a single Gaussian 
profile. However, due to the low pixel-to-pixel S/N, we cannot be 
certain of the non-Gaussianity, and hence, we decided to perform 
a single Gaussian fit. If we were to fit two Gaussian profiles, the 
resulting FWHM for both components are ~ 220 km s“^ centred 
at 93.130 and 93.223 GHz, with a total flux line similar to the one 
recovered with the one fit procedure. 


3.1.4 HLS J204314.2-214439 

In this case, we have detected four lines with S/NJS: 3.5 in the spec¬ 
trum at= 75.827, 81.178, 101.459, and 108.428 GHz. They are 
not compatible with a single redshift solution, which indicates that 
we have at least two sources inside our RSR beam. 

As a first step to associate these detections with CO transi¬ 
tions, we find that the 81.2 and 101.5 GHz lines agree very well 
with CO(4—3) and CO(5^) at « = 4.680, however there is no 
explanation for the 76 and 108 GHz lines at this redshift. The 
a = 4.680 solution has been confirmed after we reanalysed SMA 
observations, where we found the [CII] transition at t/ = 334.648 
GHz (see Fig.O. 

The Gaussian fitting to the spectral line gives a FWHM of 
432 ± 68 km s~^, which is in good agreement with the line widths 
of the CO emission lines (see TableO, and a redshift of z = 4.679 
is also in agreement with our redshift derived from the CO detec¬ 
tions. The integrated flux of the line is 19.1 ± 3.6 Jy km s“^. This 
source (hereafter HLS J204314.2-214439/S3) is the highest red¬ 
shift galaxy detected by the LMT at this time. 

Considering next the 75.827 GHz detection, we reject the 
association with CO(3-2) at « = 3.56 since the CO(4-3) line 
should appear at 101.106 GHz. The closest detected line how¬ 
ever is at 101.461 GHz, more than 1000 km/s away, and further¬ 
more we know that this line corresponds to the CO(5-4) transi¬ 
tion at z — 4.680. Hence we discard this a = 3.56 redshift so¬ 
lution. Other evidence that the 76 GHz line arises from a differ¬ 
ent source is the brightness of the line and the asymmetric line 
shape seen in Figs and |4l whilst the 81 and 101 GHz lines are 
more similar in line width and shape. We therefore identify the 76 
GHz line as CO(2-l) at a = 2.040 (hereafter HLS 1204314.2- 
214439/S 1). Moreover, the HLS team independently obtained the 
same redshift through the detection of CO(4-3) and CO(3-2) with 
IRAM30m/EMIR (Walth et al., in preparation) and Ha with Very 
Large Telescope (VLT)/Spectrograph for Integral Field Observa¬ 
tions in the Near Infrared (SINFONI, Nakajima et al., in prepara¬ 
tion). 

The remaining 108.428 GHz line does not seem to correspond 
to any transition at the redshifts of the above galaxies. By inspect¬ 
ing the available VLT/SINFONI data of RXCJ2043, Nakajima et 
al. (in preparation) found, serendipitously, a strong [OIII] emitter 


HLS J20431 4.2-21 4439/S3 



Frequency [GHz] 

Figure 5. [CII] detection of HLS J204314.2-214439/S3 using SMA obser¬ 
vations. The redshift derived from this transition, z = 4.679, is in good 
agreement with the redshift derived from the CO LMT-detected lines. This 
spectrum has been obtained after combining the signal from the three leased 
images of this galaxy. The FWHM and the redshift derived from this detec¬ 
tion is in good agreement with the values derived from the CO emission 
lines detected with the LMT. 

falling within the SCUBA-2, AzTEC, and RSR beams. The simul¬ 
taneous detection of the [OIII] doublet (A = 4959/5007 A) and 
the H/1 line confirms an optical redshift of z = 3.25. This red¬ 
shift is consistent with our 108.4 GHz line if we assume it is due 
to the CO(4-3) transition. At this redshift, we should expect to 
also detect CO(3-2) at ~ 81.3 GHz, but this is very close to our 
previously identified line CO(4-3) at 2 = 4.68 (HLS J204314.2- 
214439/S3). In order to test whether there is any evidence of CO(3- 
2) at 2 ~ 3.25, we attempt a double Gaussian fit with free cen¬ 
tral frequencies, widths and amplitudes. The best-fitting Gaussian 
profiles allowed us to recover the CO(3-2) transition at a red¬ 
shift consistent with that as predicted by the CO(4-3) transition at 
2 = 3.252 and with a similar FWHM (see Fig.|4]and Tablet^. We 
conclude that these two CO lines are most likely associated with a 
third galaxy at 2 = 3.252 (hereafter HLS J204314.2-214439/S2). 
In order to estimate the errors for the fits to the CO(3-2) transition 
we implement a bootstraping method on the noise and repeat the 
Gaussian fit 1000 times. The mean and standard deviation of the 
fits are reported in Table |3 

The detections of three components at widely different red- 
shifts (z = 2.04, 3.25, and 4.68) for this blended source confirms 
the results of theoretical works that predict that some of the blended 
galaxies are physically unassociated with typical redshift separa- 
tions Az ~ 0.9 — 1.5 dHayward et al.ll2013l : lOtwlev et aDl2015l : 
iMunoz Arancibia et al.ll2015). 

3.1.5 Fainter emission lines 

Because of the wide bandwidth of the RSR we can search for ad¬ 
ditional weaker line transitions of different species that we know 
to fall within our bandpass, given the redshifts of the galaxies, but 
which are below the higher detection threshold for blind searches. 
In this way, for HLS J204314.2-214439/S3 the atomic carbon 
[CI](®Pi Po) transition at 2 = 4.680 is expected at 86.6 GHz, 
inside the RSR spectral coverage, where we indeed have evidence 
for a detection with an integrated line S/N« 3 (see Fig. |4ll. This line 
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has a line width of ~ 420 km s“^, consistent with the line widths 
of the CO transitions, and an integrated flux of 0.9 ± 0.4 Jy km 
s“^. However, since this detection is just at a 3cr level, measure¬ 
ments derived from this line have been handled, formally, as upper 
limits. 


3.2 RSR analysis: gas properties 


Having identified the most probable line transitions and redshift of 
each source, we calculate luminosities of the galaxies u sing 
the standard relation given bv ISolomon & Vanden Boul ( l2005h : 

L'co = 3.25 X lO'ScoAV dI (1 -f z)-\ (1) 

where L'qq is the line luminosity in K km s~^ pc^, Sco AH is the 
velocity-integrated line flux in Jy km s“^, Uoha is the observed cen¬ 
tral frequency of the line in GHz, and Dl is the luminosity distance 
in Mpc. We estimate Sco AH as the integral of the Gaussian fits to 
the lines, using Monte Carlo simulations that take into account the 
errors in the Gaussian parameters (i.e. peak flux density and line 
width) to estimate the errors. 

We continue with an estimation of the H 2 mass from the mea¬ 
sured L'qq which requires two steps. First, luminosities originating 
from higher transitions (Jup > 2) must be transformed to an equiv¬ 
alent ^^CO(l-O) luminosity using a brightness ratio based on an 
excitation model. Second, once the ico(i-o) been estimated, 
the H 2 mass can be derived by the following relationship: 


Hf(H2) — a.f'co(i-o)i (2) 

where a is a conversion factor in units of Mq(K km s“^ pc^)“^. 

To transform our toI/^rln-r| ^ we have used the 

brightness ratios reported by Tcarilli & Walten ( 1201 3h for SMGs 
based on all available lite rature at the time, and a val ue of a = 0.8 
Mq(K km s“^ pc^)~^ toownes & Solomoril 1998h for the mass 
transformation. As discussed by Serieanll ( l2012h . differential mag¬ 
nification could affect the line ratios of the CO ladder, however, in 
cluster lensing it is rare to find strong gradients in the magnification 
on subarcsec ond-scales and therefo re the line ratios are usually un¬ 
affected fe.g. lDanielson et alj201lh . Furthermore, the higher-J CO 
emission has bee n shown to be mo re compact than that of lower- 
J transitions le.g. llvison etalj|201lh . which decreases even further 
the probability of differential magnification. 

The resulting ^^CO luminosities and H 2 masses are reported 
in Table |4] The mean H 2 mass of the sample is (2.0 ± 0.2) x 
10^^Mq/ fi, where ^ is the lens amplification. We caution that 
these results adopt average brightness ratios to convert to L'qq^i_q^ 
which introduce extra uncertainties, as well as our adopted value of 
a w hich is a factor of ^ 4 lower than the value used in other works 
(e.g. lThcconj_et ^|201(]|), an d couki vary with metallicity and gas 
properties ^Narayanan et alj|2012h . 

We also estimated the luminosity of the 3a detection of 
[CI](^Pi Po) in HLS J204314.2-214439/S3, finding a lumi¬ 
nosity of n^ci(i-o) ~ 5.2 ± 2.6 X 10^° K km s“^ pc^. Com¬ 
paring this value with the ^^CO transitions detected for this galaxy. 


we have line ratios of L'r 


CI(l-O) 


/^CO(4- 


3 ) 


= 0.21 ± 0.11 and 


^ci(i-o)/-^co( 5 - 4 ) “ 0.56 ± 0.29. These values are consistent 
with other measurements of SMGs and quasars and do not dif¬ 
fer significa ntly from what is f ound in l ow-redshift systems, as 
discussed by Walter et al.| ( |201 ill (see also lAlaghband-Zadeh et all 
l2013l:IFopping et alj|2014h. 

As described by Harris et al. 1 1201^ , the line luminosity-line 
width relation, Lqo(i-o) ~ G'(Aufwhm)^ (see Fig.[^, could be 



Figure 6 . ^qo( 1 —O) FWHM for lensed and unlensed 

SMGs. Red circles represent the measurements from this work and 
light red is used when we adopt the alternative redshift solution for 
HLS J045518.0+070103. Estimates derived from each detected transition 


SMGs with CO(l-O) measurements iCarilli et al 

l2010l;lHarris et al.ll2010l; 

llvison et al.ll201 it iRiechers et alj|201 it iFu et al.l 

2013|) and blue asterisks 


are SMGs with higher -J CO line measurements converted to CO(l-O) from 
The red line shows the best-fitting relation for un- 


th higher -^ 
iBothwell et alJ i2013h . 


lensed SMGs = 10^'^FWHM^ derived by IBothwell et al.l ^201 


The green 

iHarris ^ 


luares are GBT and VLA CO(l-O) detections of lensed SMGs 


squares are uB 1 and VLA 
alboilli^le et alj20l4) 


used as an estimator of the lens magnification, p, if we assume that 
the magnification only modifies the observed line luminosity (and 
not the line width). This means that ^apparent = or 


La 


apparent 

T' 


L'„ 


apparent 


a(AllFWHM)^ 


(3) 


Our measurements lie within the same region occupied by 
other lensed SMGs in Fig. which provides circumstantial ev¬ 
idence that amplification due to lensing is present. Table |4] lists 
the lens magnifications derived from this equation, adopting the 
line luminosity-line w idth relati on for unlensed SM Gs derived by 
IBothwell et ( 12013h (see also lHarrisetal.l[2012|) . We measure 
modest amplification, /r « 2 — 5, for the sample. 


3.3 SCUBA-2 + AzTEC analysis: dust properties 


The dust continuum detections allow us to estimate the dust temper¬ 
ature and IR luminosity once a spectral energy distribution (SED) 
is fitted to the data. Our determinations of spectroscopic redshifts 
(see Sectionl3.lt break the temperature-redshift degeneracy (e. g. 
iBlain et alj|2002l) . albeit there is still a degeneracy between tem¬ 
peratu re and spectral index (e.g. [Hughes et al.l l l99^ : ICasev et'd] 
|2014 . 

Assuming the dust is isothermal, the dust mass, Ma, is esti¬ 
mated from 


Mi = 


S.Dl 

(1 -I- z)KvB{v,Ti) ’ 


(4) 


where Sv is the flux density at frequency v, Ki, is the dust mass ab¬ 
sorption coefficient at v, Ti is the dust temperature, and B{vTd) 
is the Planck function at temperature Ti. The dust mass absorption 
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Table 4. Physical properties derived from the lines detected in the RSR spectra and the SCUBA-2+AzT EC photometry. Column 1: name of the source; 
column 2: observed CO luminosity; column 3: CO(l-O) luminosity con veifed using the factors repoifed by ICarilli & Walt^ j2013h ; column 4: H 2 mass 
assuming a = 0.8 Mq (K km s“^ pc^)“^ jPownes & Solorn^ll998h ; column 5: dust masses estimated from the flux densities at 850 /rm and 1.1 mm; 
column 6: total FIR luminosity; column 7: amplification factor. 


Source 

f'-^COtJup) 

(xlQlO Kkms-l pc2) 

f‘^CO{l-0) 

(xlO^'^ K km pc^) 

/rM(H2) 

(xl0l° Mq) 

pALd 

(XlO® Mq) 

A^Tfir 
(XI 012 Lq) 

Al 

HLS J045518.0+070103 

20.4 ±3.6 

31.0 ±5.5 

24.8 ±4.4 

25 ±2.0 

63l® 

r-sj 4 

HLS J173140.8+225040 

8.9 ±3.0 

13.5 ±4.5 

10.8 ± 3.6 

36 ±2.0 

54l® 

2 

HLS J204314.2-214439/S1 

32.4 ± 5.0 

38.1 ± 5.8 

30.5 ±4.7 

~ 15“ 

~ 25“ 

~ 5 

HLS J204314.2-214439/S2 

2.6 ± 1.6 

3.9 ±2.3 

3.2 ± 1.9 

~ 13“ 

~ 33“ 

~ 3*’ 


3.8 ± 1.2 

8.3 ±2.5 

6.7 ±2.0 

~ 13“ 

- 33“ 

~ 3*’ 

HLS J204314.2-214439/S3 

24.9 ±3.2 

54.2 ±6.9 

43.3 ±5.6 

~ 11“ 

~ 44“ 

~ 3*' 


9.3 ± 1.5 

23.7 ±3.9 

19.0 ±3.1 

~ 11“ 

rvj 44“ 

~ 3** 


“ The uncertainties in these values are large since the real contribution from each galaxy to the total detected (blended) flux 
is unknown. We assume equal contribution to the (sub-)mm fluxes by all components. Average of the two transitions. 


follows the same power law as the optical dep th, k, on with a n or- 
malization of Kd(850/im) = 0.15 kg“^ l lDunne et alj2003l) . 

We fit modified blackbody functions fixed at the derived red- 
shift of each galaxy, to the continuum photometry data in Table 
along with the SPIRE photometry, in order to estimate the dust 
temperature. For the HLS J204314.2-214439 blended galaxies we 
adopt a dust te mperature of T = 37 ± 8 K as that measu red in 
lensed galaxies jWeiB et al.ll2013l ; see also lHarris et alj2012l) since 
the contribution from each galaxy to the total blended flux is un¬ 
known. This value is consistent with the estimated temperatures of 
the other galaxies in our sample (e.g. 38 ± 3 and 35 ± 3 K, for HLS 
J045518.0-I-070103 and HLS J173140.8-I-225040, respectively). 

Using this temperatures and a fixed emissivity index of /3 = 
1.7 (empirical fits to the observed long wav el ength SEDs suggest 
p = 1.5 — 2; e.g. iDunne & EalesI 1 1200 ih ; IChapin et alj 1 20091 ; 
iMagnelli etZI|2012|) . we have derived the dust mass for each 
galaxy from the flux density at both 850 /rm and 1.1 mm, and 
extrapolating k to the observed rest-frame for each galaxy. These 
values are reported in Table |4] and have a mean dust mass of 
(2.0 ± 0.2) X 10® Mq/zt. If we extrapolate the observed spectral 
index to use Kd(850/im) directly, the mean dust mass increases 
by a factor of ~ 2. These calculations do no t include uncertain - 
ties in k, which could be at least a factor of 3 tounne et al.ll2003h . 
For HLS J204314.2-214439/S1, HLS 1204314.2-214439/S2, and 
HLS J204314.2-214439/S3 we have assumed they each contribute 
a third of the total (blended) flux density, as a first approximation. 
A full multi-wavelength analysis of this triple system, including 
interferometric observations, will be presented by Walth et al. (in 
preparation). 

Finally, the estimated FIR luminosites (Lfir; 42.5 — 
122.5/rm), which arise from the SED fitting, are reported in Table 
|4]and discuss below. 


4 DISCUSSION 

4.1 The V!co "L fir relation 

The SFE, with which the molecular gas is being turned into stars, 
can be inferred from our observations. This quantity is often es¬ 
timated by the ratio SFR/M(H 2 ) - the inverse of the gas deple¬ 
tion time - however, a more straightforward estimate of the SFE is 
the continuum-to-line luminosity ratio, I/fir/Lco(i-o) > because 


it does not depend on the CO-H 2 conversion factor, a. This de¬ 
scribes, in observational terms, the relationship between the lumi¬ 
nosity due to star formation and the gas content. 

In Fig. [ 7 ] we plot our galaxies with CO detections, after 
correcting for magnification, on to the L'co — Lfir plane, in¬ 
cluding obs ervations of local (ultra-) luminous infrar ed galaxies 
(U)LIRGs I Sanders^t^ j^9ll ; ISolomon et al.l 1 19971) and other 
SMGs teothwell et al.ll20 13l). We al so show the three power-law 
fits derived bv iBothwell et al.l ( 1201 3h to the local (U)LIRGs alone 
(with a slope of 0.79 ± 0.08), the SMGs alone (0.93 ± 0.14), and 
the combined sample (0.83 ± 0.09). Our results are consistent with 
other SMG observations, but they are also consistent with the fit to 
the combi ned sample and the (U)LIRGs alone. Furthermore, as dis¬ 
cussed bv iBothwell et al.L our analysis requires extrapolating from 
high-Jup ^®CO transitions to ^^CO(l-O), which introduces extra 
uncertainties without brightness temperature ratio measurements 
for individual sources. Although we have used the estimated am¬ 
plification factors, the ratio between CO and FIR luminosity is in¬ 
dependent of magnification, assuming no differential amplification, 
which is usually the case in the cluster lensing. 


4.2 The gas-to-dust ratio 


Constraining the molecular gas properties of a large sample of non- 
gravitationally amplified SMGs with spectroscopic observations re¬ 
quire s a significant investment of time at present teothwell et al.l 
l2013h . An alternative approach to estimate the molecular gas mass 
is to use the continuum emission of the dust in order to estimate the 
dust mass, and then use an appropriate gas-to-dust ratio, 5gdr, to 
derive the total mass of the molecular ISM. However, this gas-to- 
dust ratio is only well constrained in Mil ky Way molecular clouds 
and local galaxies (e. g. Draine et al .l2007h . and only in a few SMGs 
(e.g. lSwinbank et S]|2013h . 

Combining our estimation of H 2 gas mass with the dust mass 
resulting from the SCUBA-2 and AzTEC flux densities, we esti¬ 
mate the gas-to-dust ratio for these high-redshift galaxies, which 
is independent of the lensing amplification, assuming that both 
gas and dust are amplified by the same factor, which should be 
applicable to these cluster-lensed sources. Using our independent 
estimation of these quantities we have found a weighted mean 
gas-to-dust ratio of 5 g dr ~ 55 for the two galaxies without 
direct evidence of blending (HLS J045518.0-I-070103 and HLS 
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Figure 7. Lens-coiTected CO luminosity versus FIR luminosity {Lqq — 
Lfir) for our sample, in red circles, and in light red if we adopt the al¬ 
ternative CO solution for HLS J045518.0+070103. For galaxies with more 
than one detected transition we have plotted the m ean For compari¬ 

son, other SMGs are plotted iBothwel l et alJl oT^ alon g with local LIRGs 
iSanders et al.|l9 9l[) and U LIRGs isolomon et al.ll997h . Best-fitting slopes 
derived bv lBothw*^ et alj l2013l) to SMGs, local (U)LlRGs, and all three 
combined samples are overplotted. They have slopes of 0.93 ± 0.14, 
0.79 ± 0.08, and 0.83 ± 0.09, respectively. 


J173140.8+225040). If we also consider the blended galaxies, as¬ 
suming they all contribute equally to the total integrated flux with 
intrinsic flux density errors of 50 per cent, we find a weighted 
mean (5gdr ~ 75. The scatter in these measurements is how¬ 
ever high ((5gdr = 20 — 390), especially for the blended sources. 
Our estimated mean values are consistent with tho se of non lensed 
SMGs, 5 gdr = 90 ± 25 dSwinbank et al.ll2013h . and also with 
the le nsed 2 = 6.3 HFLS3 galaxy, 5gdr ~ 80 dRiechers et al.l 
l2013h . These ratios are s maller than local values, as also dis¬ 
cussed in the literature (e.g .jKovacs et alj200^ : ISantini et al.l2O10l : 
IPessauees-Zavadskv et al.ll2015h . which supports the ‘dust rich¬ 

ness’ interpretation of SMGs. For reference, the Milky Way has 
a gas-to-dust rati o of ~ 135 — 185 (Draine et al. 2007), and local 
LIRGs 120±28 dWilson et alj200^ . We emphasize, however, that 
the scatter in our measurements, especially for the blended sources, 
is high, and some of our estimated individual ratios for the blends 
overlap with the local values. 


5 SUMMARY 


We present AzTEC 1.1 mm continuum observations and RSR spec¬ 
tra taken with the LMT, as part of the LMT Early Science cam¬ 
paign, of a sample of four cluster-lensed SMGs. 

We surprisingly find that one of the four targets studied is a 
blend of three galaxies at different redshifts. This confirms that the 
multiplicity that has been found by high-resolution interferomet¬ 
ric observations in norma l SMGs dWang et al.ll201lLISmolcic et al.l 
I2OI2I : iHodge et alj l2013h is also present in lensed galaxies, al¬ 
though our search focuses on spectroscopically distinct compo¬ 
nents. This also confirms the theoretical prediction s of physically 
unassociated components in blended galaxies le.g.lHavward et al.l 

l2013LICowley et^l2015l : lMunoz Arancibia et al.l2015l) . The sam¬ 


ple, though, is too small to derive firm conclusions on the statistics 
of spectroscopic multiplicity. 

Of the five reported redshifts, just one was obtained with the 
robust blind detections of two lines within a single RSR spectrum. 
Three other redshifts were confirmed thanks to the detection of dif¬ 
ferent transitions with other telescopes, and the remaining one is a 
tentative solution. 

The estimated mean H 2 gas mass estimated from the CO lines 
detected is (2.0 ± 0.2) x 10^^ Mq// i and the mean dust mass esti¬ 
mated from the continuum data is (2.0 ± 0.2) x 10 ®Mq/^, where 
is the lensing factor. These results are self-consistent using either 
the 850 /rm or 1.1 mm data. Using the lin e width-luminosity rela¬ 
tion for SMGs from fBothwell et all (1201 3l) we estimate ^ 2 — 5. 

We infer from our independent estimations of gas and dust 
masses a weighted mean gas-to-dust ratio (5gdr « 55 — 75, which, 
given the little likelihood of differential amplification at mid-J tran¬ 
sitions, should be independent of the amplification factor. These 
values are consistent with other measurements of SMGs and lower 
than the estimated ratios for local galaxies. 

Einally, we estimate the FIR luminosity using our photometric 
data and with this the SFE of our galaxies through the L'qq — Lpm 
relationship. Our bright galaxies follow the correlation found for 
local (U)LIRGs and other SMGs, with a close to linear slope. This 
suggests that, at the largest luminosities sampled by our systems, 
the SFE is maintained arguing for a universal value from the lowest 
to the highest redshifts and from the smallest to the largest mass 
reservoirs of molecular gas. 

The successful blind redshift search in these bright SMGs, per¬ 
formed during this LMT Early Science operation period with the 
current 32-m diameter illuminated surface, highlights the compet¬ 
itiveness of this new facility. The sensitivity that will be achieved 
with the expanded primary surface of the full 50-m diameter LMT 
will place this telescope among the m ost powerful facilit ies to 
perform spectroscopic CO surveys (see iFraver et al.l 1201 ll for a 
comparison of different telescopes), allowing us to estimate red¬ 
shifts for significant samples of highly obscured dusty star-forming 
galaxies which are not measurable with even the largest optical and 
near-IR telescopes. 
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